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The problem of the nuclear spin-lattice relaxation in La2Cu04 is revisited in con- 
nection with the recent measurements of the NQR relaxation rate for temperatures 
up to 900K [T. Imai et al, Phys. Rev. Lett., in press]. We use an approach based on 
the exact diagonalization for the Heisenberg model to calculate the short wavelength 
contribution to the relaxation rate in the high temperature region, T •> J/2. It is 
shown that the spin diffusion accounts for approximately 10% of the total relaxation 
rate at 900K and would become dominant for T>J. The calculated 1/Ti is in good 
agreement with the experiment both in terms of the absolute value and temperature 
dependence. 
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The fact that the spin dynamics of the parent insulating compound La2Cu04 is described 
by the 5 = 1/2 Heisenberg model with J~ 1500K is now very well established (for reviews, see 
Recently, T. Imai et al. have measured the copper nuclear spin-lattice relaxation 
rate, \jT\ in the undoped and Sr-doped L^CuCU for temperatures up to 900K. They find a 
plateau in 1/Ti as a function of temperature for 700<T<900K. In this temperature region, 
the relaxation rate is insensitive to doping, a result which suggests that at high temperatures 
the dominant relaxation mechanism is the same in both metallic and insulating samples |§ . 

As it is known, for localized spins, the relaxation rate is determined by the so-called "ex- 
change narrowing" mechanism ||]. The "exchange narrowing" here refers to the relaxation 
process governed by the spin-spin exchange interaction. An approach based on the Gaussian 
approximation for the dynamic structure factor has been developed in Ref. |3| in order to 
calculate the relaxation rate for T^> J. In Ref. [||,[| , this approach has been combined with 
high temperature expansion method and thus extended to finite temperatures of the order 
of J. For temperatures larger than J, 1/Ti has been shown to increase as the temperature 
increases. On the other hand, in the low temperature limit the dominant contribution to the 
copper relaxation rate is due to critical fluctuations around q=(ir/a,7r/a), and it increases 
exponentially as the temperature decreases, 1/Ti oc T 3 / 2 exp(27rp s /T) |7j]. For T< J, the 
spin stiffness is p s ~ 0.18J. The interpolation from low to high temperatures shows that 
1/Ti as a function of temperature has a minimum. In Ref. |7j], its position has been predicted 
at T ~ 700K, a result which seems to be in contradiction with the experimental data of 
Ref. §. 

Therefore, the purpose of this work is to understand whether or not this experimental 
result can be quantitatively understood in the framework of the nearest-neighbor Heisenberg 
model. The analysis of the NMR data in La2_ a; Sr a ,Cu04 has lead to the conclusion that the 
hyperfine constants in this material approximately coincide with those of YB^CusOz ||. 
We take advantage of this and use the values of the hyperfine couplings obtained in Ref. 
H for the yttrium-based compounds. Along with the use of J ~ 1500K for the exchange 
constant, this eliminates all adjustable parameters in our calculation. 
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The copper spin-lattice relaxation rate measured in the NQR experiment is: 

1 = ^ linj *M, (1) 

where 

(for simplicity, we use the units where ks — h — 1). In the NQR experiment, the hyperfme 
formfactor A(q) is given by |§: 

A(q) = A xy + 2B cos(q x a) + 2B cos(q y a), (3) 

where A xy and S are the in-plane local and isotropic transferred hyperfme couplings, re- 
spectively. In what follows, we use A xy /B = 0.84, B = 40.8KOe//i B [§. 

The relaxation rate of the Heisenberg antiferromagnet has been discussed in several 
publications. However, the low temperature calculation based on the dynamical scaling 
theory |7|] is not valid for T~ J/2 = 750K, where also the contribution from wave vectors 
other than q— (n/a, n/a) becomes important. On the other hand, it has been mentioned in 
Ref. H that the high temperature expansion results based on the Gaussian approximation 
do not show the low temperature increase of 1/Ti, apparently because of the particular 
functional form assumed in this calculation for the dynamical structure factor. Using large 
N expansion technique, it has been shown in Ref. JTO] that 1/Ti is nearly temperature 
independent for T ~ J/2. Unlike our calculation, this approach does not start from the 
S = 1/2 lattice model and so the absolute value of 1/Ti is evaluated in Ref. using the 
low-temperature fit of the same data and not the hyperfme couplings A xy , B. 

For T <^ J, the spin diffusion (q — > 0) contribution to the relaxation rate, (l/Ti)^//, 
is negligible because the spin diffusion constant, D, is exponentially large pynfl. However, 



D rapidly decreases as the temperature increases, that is, the q — > component may be 
important for higher temperatures. In a pure two-dimensional model, the conservation of 
spin leads to the divergence of (l/Ti)^/; that is, the relaxation would be faster than expo- 
nential. However, in a real system (1/Ti)^// remains finite and its magnitude is determined 



by the length scale L s , set either by spin-nonconserving forces or three-dimensional effects. 
Since in any cluster calculation (exact diagonalization or Monte-Carlo) the cutoff is set by 
the lattice size, we have taken into account the q— >0 contribution separately. 

Our approach for the calculation of the short wavelength contribution to the relaxation 
rate is based on the exact diagonalization of the Hamiltonian for the 4x4 cluster. Since the 
nuclear spin-lattice relaxation rate is determined by short-range spin correlations, our results 
are relevant to the real system as long as the correlation length is not large compared to the 
cluster size. The spectral representation for x'hf can be written in terms of the Hamiltonian 
eigensystem as follows: 

Xfc/H = £ 
7 

ab 



, , ~ 7 E bM-E a /T) - exp(-E b /T)] 



x b{E a -E h ^)W^\^)MS z q \h)\\ (4) 

q^O 

where E a ^ are the eigenvalues of the Hamiltonian and Z = J2a ex P{~E a /T) is the partition 
function. In the thermodynamic limit (AT— >-oo), x'hf * s a continuous function of frequency, 
while for finite size it is a superposition of delta functions. 

For a finite cluster, the limit u — > in Eq. ([I]) is not defined, but we argue that the ther- 
modynamic x'hf(u) can be calculated using the following procedure. Consider the auxiliary 
function given, for a cluster of size N, by 

.2, .2 r r ,A _ 1 r a, rtfie) 



g fi B I N (u) = - de —± — . (5) 

From this equation, Xhfi^)/^ = Q 2 l^B^dl ^ / dui) . For finite cluster, In(uj) can easily be 
calculated from the eigenstates of the Hamiltonian: 

I N (w) = J^Iab [B{E a -E h +u)-B{E a -E h -u)\ , (6) 

ab 

where 9{x) is the Heaviside function and 

7i exp(-E a /T) - exp(-E b /T) 



1 ab 



2Z El — E a 

1 

N 



x^E^(q)l(«l^l^)l 2 - (7) 



(for a = b we take the limit E a ^Eb). The auxiliary function In(u) is quite smooth as long 
as the temperature is not much smaller than the gap between the ground state and the rest 
of the spectrum, which for the 16-site cluster is of order J/2. For temperatures T> 1.5—2 J, 
we find no appreciable size dependence: iio — he- In the study of static properties of the 
Heisenberg model [|12|| , no discrepancy was found between the 4x4 cluster and Monte-Carlo 
results for larger systems at T > J. Both the discrepancy and the errorbars in the fitting 
of ii6 by a smooth function increase up to approximately 10% for T~ J/2. Thus, we will 
assume that our calculation of the short wavelength contribution to 1/Ti has 10% accuracy. 

Now we turn to the calculation of the q — > contribution to 1/Ti. For T" 1 < q < 
max(£,a) _1 and lot <C 1, the dynamical spin susceptibility x(q, cu) has the following form: 

where D is the diffusion constant and r a characteristic relaxation time. Substituting this 
expression into Eq.(|I]), we obtain 

1 \ T X oa 2 A 2 (q = 0) 



r I 2 2 n log 7 ' ( 9 ) 

-TlSdiff ^9 ^B D L f.s. 

where we take J/. s . > £ to be equal to the size of our cluster. For T ^> J, the diffusion 
constant is D ~ 0.43Ja 2 fL3|, so that (l/Ti) diff ~ 7400 x \og(L s /L f , s ) sec" 1 is at least 
several times larger than the measured rate at the maximal accessible temperature 900K || . 
This contribution is larger than the calculated short wavelength contribution at the same 
temperature. Therefore, the relaxation rate of the 2D Heisenberg model for T> J is a poorly 
defined quantity, since it strongly depends on the way the cutoff is taken into account. In 
this temperature region, an accurate calculation of 1/Ti would have to involve the actual 
mechanism destroying the diffusion. However, since temperatures larger than the exchange 
constant are not experimentally accessible in La2Cu04, we will examine now whether or not 
the spin diffusion substantially contributes to the relaxation rate at 600 — 900K. 

In order to address this issue, we have to determine both the diffusion constant D and 



the length scale L s for T~ J/2. The diffusion constant is estimated as |14|,|15 



D = 7 limg- 2 (^)f <<4>-V2, (10) 
where (k>2 n ) g are the frequency moments of the dynamical response function, 



and 7 is a numerical factor which depends on the assumed short time relaxational behavior 
P13| ]. Taking 7 = ps 1.25 [|15| yields D T=0 o = OAOJa 2 , which is quite close to the value 



0. 43 J a [lTSJ obtained through an evaluation of the memory function. In Ref. |]17||, general 



expressions for the series in (3 = J/T for (cf 2Tl ) g have been derived. Using these results, we 
calculate first two terms of the high temperature expansion for the diffusion constant: 

D _ v^F + 2Vp? + Q(/?2) ^ Q 4Q + QA2j3 + 0(/?2) _ (12) 



Ja 2 2^ 40^ 

Two leading terms in the high temperature expansion series are not sufficient for the accurate 
estimate of D at T ~ J/2. However, we know that the diffusion constant should scale 
approximately as D oc £ in the quantum critical region, p s <T <J. Using Monte-Carlo data 
of Ref. [JTSJI for the correlation length, £(T), we estimate D ~ 3Ja 2 at T = 900K, which when 
substituted into Eq.(|) gives {1/T 1 ) difl ~ (200-300) x log(L s /L/. s .) sec" 1 . 

Now we turn to the evaluation of the logarithm in Eq.(|S]). Since the hyperfine splitting ~ 
1.5 • 10~ 7 eV is very small and above the tetragonal-to-orthorhombic transition temperature, 
Tt-o — 525K, the Dzyaloshinskii-Moriya interaction vanishes, the cutoff is determined either 
by the three dimensional effects or by the nonconservation of spin. Consider first the cutoff 
due to the three-dimensional effects, L^ , which is set by the interplanar diffusion constant, 
-D_i_. For estimation purposes, we express D in terms of the characteristic damping of spin 
waves for small wave vectors, 7 03 • With omission of all factors of the order of unity, we 



get D±/D\\ ~ J'/ J, which yields L^ D ~ 300a, a quite large value. 

Given the size of Lg D , we consider an alternative physical origin for the cutoff, the 
presence of weak disorder in Cu0 2 planes. For temperatures above 700K, the oxygen content 
changes after the heating cycle by approximately 0.004 per unit cell 0; that is, the average 
distance between nonstoechiometric oxygen atoms, which we identify with L s , is 10 — 20a. 



Although the value of L s cannot be determined quite accurately, the q — > contribution 
to 1/Ti depends on L s only weakly. In what follows, we plot the results for L s = 10a 
and L s = 20a. Substituting the above values of L s into Eq.(||), one obtains that the spin 
diffusion contribution accounts for approximately 10% of the total spin lattice relaxation 
rate for T = 900K, but rapidly decreases as the temperature decreases. This explains why 
the tetragonal-to-orthorhombic transition at Tt-o — 525K does not have any observable 
effect on the spin lattice relaxation although it affects L s . 

The total relaxation rate for L s = 10 and L s = 20 and the short wavelength contribution 
alone are plotted in Fig.|l] together with the experimental result of Ref . || . The theoretical 
result is in 15% agreement with the experiment (for L s = 10). The agreement can be 
improved by either taking smaller L s , or changing the hyperfine constants by 7% (actually, 
A xy ,B are known only with 5-10% accuracy ||). It is important to emphasize, however, 
that the ambiguity in definition of 1/Ti as a function of the cutoff L s exceeds our estimate 
of the systematic error of the finite cluster calculation; moreover, the cutoff itself can not 
be determined quite accurately. The spin diffusion (q — > 0) contribution rapidly increases 
as the temperature increases (Fig.|I|, inset) and becomes dominant for T > 1.5J, as it is 
shown on FigfJ. Although this temperature range is beyond the limit of chemical stability 
for La2Cu04, it may be of interest for other materials described by the Heisenberg model 
but with smaller J, such as Cu(HC0 2 ) 2 -4H 2 and Cu(pyz) 2 (C10 4 ) 2 @- 

To summarize, we have calculated the copper spin-lattice relaxation rate for La 2 Cu04 
without introducing any adjustable parameters. The spin diffusion (q — > 0) contribution is 
shown to account for 10% of the relaxation for the maximal temperature achieved in the 
experiment, 900K, although it would become dominant for larger temperatures, thereby 
explaining the discrepancy between different calculations of the relaxation rate at high 
temperatures. The measured 1/Ti turned out to be quantitatively consistent with the 
nearest-neighbor Heisenberg model description of the spin dynamics in La 2 Cu04. 
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FIGURES 

FIG. 1. The calculated 1/Ti as a function of temperature without (solid line) and with the 
spin diffusion contribution for L s = 10 (dashed line), and for L s = 20 (dotted line). Dots are the 
experimental result of Ref. ||. The errorbars due to the calculation inaccuracy (less than 10%, 
not shown) are smaller than the ambiguity in the definition of 1/Ti related to q — >0 cutoff. Inset: 
the spin diffusion (q— >0) contribution to the relaxation rate. 

FIG. 2. Same as Fig.|l|, but in the temperature range J/2 < T < 3J for a hypothetical 
heat-resistant sample. 
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